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Resume 

On s'interesse dans cette note a I'equation de Fisher-KPP dans I'espace entier, 
ou le laplacien est remplace par le generateur d'un semi-groupe de Feller a noyau 
lentement decroissant, un exemple important etant le laplacien fractionnaire. A 
la difference de I'equation de Fisher-KPP classique, oii I'etat stable envahit I'etat 
instable a une vitesse constante en temps, nous montrons que la vitesse d'invasion est 
exponentielle en temps. Ces resultats apportent une justification mathematiquement 
rigoureuse a de nombreuses heuristiques sur ce modele. 

Abstract 

Front propagation in Fisher-KPP equations with fractional diffusion. We 

study in this note the Fisher-KPP equation where the Laplacian is replaced by the 
generator of a Feller semigroup with slowly decaying kernel, an important example 
being the fractional Laplacian. Contrary to what happens in the standard Laplacian 
case, where the stable state invades the unstable one at constant speed, we prove 
here that invasion holds at an exponential in time velocity. These results provide a 
mathematically rigorous justification of numerous heuristics about this model. 
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Abridged English version 



Consider the reaction-diffusion equation 

Ut + Au^f{u){t>0,xeR^), u{0,x) e[0,l]. (1) 



where A is the generator of a Feller semigroup (see [3]). That is, the solution 
of Vf + Av — 0, ^(0, x) — Vo{x) is given by 

e-*^Vo = P{t, ■) * Vq, 

where p is a positive kernel satisfying assumptions (i)-(iv) below. In the whole 
note we will assume 

^ <p{t:x)<^__ , ^1 , 0<a<l,(2) 



where S is a constant. An important example for A is the fractional Laplacian 

/ A / \ f u(x) — u(y) , 



where the above expression has to be understood in the sense of principal 
values. The function / is concave, and satisfies /(O) = /(I) = and /'(I) < 
< /'(O). When A — —A (i.e. a — 1), equation (1) is known as the Fisher- 
KPP equation [7]. 

We wish to understand the large-time behaviour of (1). In view of the classical 

reference [1] - which treats the case a = 1 - we expect the stable state 1 to 
invade the unstable state 0, and we wish to know at which rate. In [1], it is 
linear in time. Numerous heuristics from physics papers - see for instance [8] 

- suggest that, in the case a e (0, 1), the invasion rate will be exponential in 
time. Our results provide a mathematically rigorous proof of this. 

The first result concerns compactly supported initial data. 

/'(O) 

Theorem 0.1 Let a G (0, 1) and set c* = — - — — . Assume u(0, ■) ^ to be 

^ ^ N + 2a V w ^ 

compactly supported. Then: 

- for all c < c*, we have lim inf u(t, x) — 1, 

t->+oo |a;|<e'=* 

- for all c > c*; we have lim sup u{t, x) — 0. 
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The following result deals with nondecreasing initial data. 

f'(0) 

Theorem 0.2 Let a e (0, 1) andN =1. Set c„ = Assume u{0, •) ^ 

2a 

to be nondecreasing and such that M_ fl supp(m(0, •)) is compact. Then: 

- for all c < c**, we have lim inf u(t,x) = 1, 

- for all c > c**, we have lim sup u{t,x) — 0. 

The proofs of these theorems also give that (1) cannot have travelling wave 
solutions if a < 1, in sharp contrast with the case a — 1. 

In the case a — 1/2 and N — 1, we have sharper asymptotics on the dynam- 
ics of the level sets of solutions with compactly supported initial data as in 
Theorem 0.1. 

Theorem 0.3 Let a — N — 1, and f{u) = u{l — u). Assume m(0, •) ^ 
to he compactly supported. For X e (0, 1), set 

x'^{t) — mf{x : u{t,x)—X} and xj^{t) — sup{x : u{t,x) — X}. 

Then, there exists a constant C\> 1 such that, for t large enough, 

-CAe*/2 < < --^e*/2 and -^e*/^ < < Cab^/^ 

All these results will be established and developed in [4]. 
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1 Introduction, motivation 



Soit A le generateur d'un semi-groupe de Feller (voir [3]), i.e. un operateur sur 
I'espace des fonctions continues bornees sur tel que le probleme d'evolution 

vt + Av = {t> 0, X eR^), v{0,x)=vo{x) (3) 

admette une unique solution v{t,x) pouvant etre representee sous la forme 

v{t,x) = j p{t,x - y)vQ{y) dy, 

Oil le noyau p verifie les proprietes suivantes : 

(i) . p{t, a;) > et p{t, •) est de masse unite, i.e. j p{t, x) dx = 1; 

(ii) . p reguliere sur ]R_|_ x ; 



(iii). p{t, •) * p{s, •) = p{t + s, •) pour tout couple (s, t) e 



(iv). il existe S > 1 tel que, pour tons (t, x) e 1R+ x on ait : 

^ <p(i,2:)<^^- , , 0<a<l.(4) 



Notons que si A est le generateur d'un processus de Markov stable, i.e. il cxiste 
une fonction a(t) telle que p(t, x) = a(t)~^p(l, x/a(t)), alors il existe a G (0, 1] 
tel que a{t) = - comme dans I'expression (4). Un exemple important d'un 
tel operateur est le laplacien fractionnaire (—A)", avec < a < 1 : si est 
assez reguliere et a croissance assez lente a I'infini, 

i-i^Tuix) = c^,„ / dy, 0<a<l, 

oil I'expression ci-dessus est a prendre au sens des valeurs principales, et oii 
la constante catq, est ajustec pour que le symbole de (—A)" soit Voir 
par exemple [5]. Le cas a = 1/2 est un cas important oii p{t,x) admet une 
expression explicite : 

p{t,x) 



{t^ + |x|2)(^+l)/2 t^{l + |i-lx|2)(^+l)/2 ' 

Oil est choisie pour assurer la propriete (i). 



4 



Soit / unc fonction dc classe sur [0, 1], concave, telle que /(O) = /(I) = 
et /'(I) < < /'(O) - par exemple f{u) = u{l — u). On s'interesse au 
comportement en grand temps des solutions u{t,x) du probleme d'evolution 

ut + Au^f{u){t>0,xeR^), u{0,x) e[0,l\. (5) 

Quand A = —A (i.e. a = 1), (5) est connue sous le nom d'equation de Fisher- 
KPP [7]. Le probleme de Cauchy pour (5) est, de par les proprietes (i)-(iii) 
de A, trivialement bien pose sur M+ x , et verifie de plus le principe de 
comparaison suivant : pour deux solutions u et v dc (5), si w(0, •) < ^'(0, ■) 
alors u < V sur IR+ x R^. On s'attend a ce que I'etat stable 1 envahisse I'etat 
instable 0, et on souhaite savoir a quelle vitesse. 

Ce type d'equation de reaction-diffusion apparait dans de nombreux modeles 
de physique (turbulence, plasmas chauds, modeles de flammes) lorsque les 
phenomenes diffusifs ne sont pas correctement dccrits par des processus gaus- 
siens. Voir par exemple [8] pour une description de quelques-uns de ces modeles. 
Bien qu'heuristiques, ceux-ci peuvent apporter des informations qualitatives. 
L'equation (5) apparait aussi dans une classe de modeles de dynamique des 
populations, et pent s'obtenir, dans un certain regime spatio-temporel, comme 
I'asymptotique d'un modele integro-differentiel - voir [2]. L'equation (5) avec 
A = —A (i.e. a = 1) pent s'obtenir comme asymptotique du meme modele, 
dans un regime spatio-temporel different - voir [7] . II est done legitime de cher- 
cher a comprendre les proprietes qualitatives de (5). De fait, de nombreuses 
heuristiques physiques existent, voir la encore [8]. Le but de ce travail est 
d'apporter une justification rigoureuse sur le plan mathematique. 



2 Resultats principaux 

Lorsque a — 1 (equation de Fisher-KPP) on rappelle que l'equation (5) devient 
ut- Au^ f{u) {t>Q, X eR^), m(0,x) e [0,1]. (6) 

Un rcsultat tres general sur (6) est du a Aronson- Weinberger [1] ; il decrit le 
devenir d'une perturbation a support compact. 

Theoreme 2.1 ([1]) Soit u solution de (6) et supposons que u{0, •) ^ est 
a support compact. Conside-rons la quantite c*^i = 2^/'(0). Alors : 

- pour tout c < c* 1, on a lim inf u(t, x) — 1, 

t-*+oo\x\<ct ^ ^ 
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- pour tout c > c* 1, on a lim sup u{t, x) — 0. 

II se trouve que (6) admet des solutions d'ondes progressives monodimension- 
nelles connectant a 1, i.e. des solutions de la forme (f){x ■ e + ct), avec 

-0" + c0' = 0(-oo) = 0, <^(+oo) = 1. 

Le reel c*^i est justement la plus petite vitesse d'onde possible, et Komogorov, 
Petrovskii et Piskunov [7] ont demontre que la solution de (6) avec = 1, 
issue de la fonction de Heaviside H{x), converge en grand temps vers un profil 
d'onde de vitesse c*_i. 

Nos resultats sont les suivants. 

/'(O) 

Theoreme 2.2 Soit a G (0, 1), et considerons la quantite c* = . Soit 

^ ^ ^ N + 2a 

u solution de (5) et supposons que u{0, •) ^ est d support compact. Alors : 

- pour tout c < c*, on a lim inf u(t,x) — 1, 

- pour tout c > c*, on a lim sup u{t, x) — 0. 

Ce type d'invasion exponentielle a ete remarque par Hamel et Roques [6] 
dans le cadre de I'equation (6) (i.e. quand A est le laplacien standard) pour 
des donnees initiales plus grandes que une fonction positive a decroissance 
algebrique vers a I'infini. 

Notre resultat suivant concerne les solutions avec donnee initiale croissante. 

/'(O) 

Theoreme 2.3 Soit a G (0, 1), N = 1, et considerons la quantitec^^ = . 

2a 

Soit u solution de (5) et supposons que m(0, ■) ^ est croissante et telle que 
]R_ n supp(m(0, •)) soit compact. Alors : 

- pour tout c < c**, on a lim inf u(t,x) — 1, 

- pour tout c > c**; on a lim sup u{t,x) = 0. 

On montre en particulier que (5) ne pent pas admettre d'onde progressive si 
a < 1, en contraste avec le cas a — 1. 

On remarquera une difference supplement aire avec le cas « = 1 : le theoreme 
de KPP pour a = 1 assure que la vitesse d'invasion de I'etat stable dans le cas 
d'une donnee initiale croissante est la meme que pour une donnee a support 
compact. Dans le cas a < 1, on a deux vitesses differentes. 
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Le dernier resultat concerne une asymptotique plus precise de la dynamique 
des lignes de niveau. 

Theoreme 2.4 Soient a = -, N = 1 et f{u) = u{l — u). Soit u solution de 
(5) et supposons que m(0, •) ^ est a support compact. Pour \ e (0, 1), soient 

x'^{t) — vai{x : u{t,x) — X\ et x'^{t) — svlp{x : u{t,x) — X\. 

Alors, il existe une constante Cx> 1 tel que, pour tout t > assez grand, 

-Cxe"'<x-^{t)<-^e'l' et ^e^'^ < xlit) < C^e^'^ 

Ce resultat permet de poursuivre I'etude de la dynamique complete de (5). II 
admet sans doute une generalisation pour tout a e (0, 1), que nous ne savons 
pas effectuer a I'heure actuelle. 



3 Elements de preuves 

Dans le cas a — une fagon de faire apparaitre la quantite c*,i est la suivante : 
puisque / est concave, une sur-solution de (6) est la solution u{t, x) de 

Ut-Au = f{0)u, u{0,x) =u{0,x). 

Examinant le cas particulier u{0, ■) — Sq - la masse de Dirac en - on ob- 

tient u{t,x) = (47rt)^e-'^'(°)*-^, et u ^ X si \x\ = 2^/'(0)t + o{t). Ces 
considerations tres simples permettent aussi de comprendre pourquoi une 
decroissance algebrique de la donnee initiale - comme dans [6] - conduit aussi 
a une invasion exponentielle en temps. 

Utilisons cctte heuristique - deja effectuee dans [8], par exemple - dans le cas 
a < 1. Cette fois-ci, la solution u{t,x) de 

ut + Au = f{0)u, u{0, ■) = 6o 

est 

u^ef'mp^t,x). 
L'estimation (4) indique que TZ = A si 

\x\ = 0(ti^)e-*, c. = -f^. (7) 
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La propriete de sur-solution donne une borne supcricure sur la vitesse des 
ensembles de niveaux et, partant, la deuxieme partic du Thcorcmc 2.2. L'ob- 
tention - plus difficile - de I'autre partie (i.e. la convergence vers 1) se fait 
en deux temps. D'abord, on montre I'existence de £ G (0, 1) tel que, si cr < 
f'{0)/{N + 2a), alors 

liminf inf u(t,x)>e. (8) 



Ceci est accompli en construisant de solutions de I'equation Vt+Av — S~^f{S)v 
qui prennent des valeurs en (0, 5) - qui sont en consequence des sous-solutions 
de (5). Pour ceci, on resout I'equation pour v avec une donnee initiale prenant 
des valeurs dans (0, e) (avec e < 6) jusqu'au temps T oii la solution prend 
la valeur S, puis en coupant v{T, •) par e, et en utilisant de fagon recursive 
cet argument. La convergence vers 1 est prouve en utilisant (8) et une sous- 
solution liee a la solution de I'equation homogene (3) avec donnee initiale \x\'^, 
7 e (0,2a). 

La demonstration du Theoreme 2.4 sur les ensembles de niveaux des solutions 
repose sur des arguments similaires combinees avec I'utilisation des sous et 
sur-solutions plus precises - dans ce cas Q; = l/2etA^ = l-dela forme 

pour certaines constantes a et fonctions b = b{t) du temps. Notons que ce 
Theoreme 2.4 etablit que la position du front donnee par I'argument heuris- 
tique (7), i.e. \x\ — 0(i^/^)e*/^, n'est pas correct, et a la place est donnee par 
\x\ = 0(l)e*/2. 

Tons ces points seront developpes dans [4]. 
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